(12) NACH DEM VERTRAG UBER DIE INTERNATIONALE ZUSAMMEN ARBEIT AUF DEM GEBIET DES 
PATENTWESENS (PCT) VEROFFENTLICHTE INTERNATIONALE ANMELDUNG 



(19) Weltorganisation fiir geistiges Eigentum 

Internationales Biiro 

(43) Internationales Veroffentlichungsdatum 
28. Juli 2005 (28.07.2005) 




PCT 



IIH 

(10) Internationale Veroffentlichungsnummer 

WO 2005/068953 Al 



(51) Internationale Patentklassifikation 7 : 



G01J 3/46 



(21) Internationales Aktenzeichen: PCT/EP2005/0003 1 9 

(22) Internationales Anmeldedatum: 

14. Januar2005 (14.01.2005) 



(25) Einreichungssprache: 

(26) Veroffentlichungssprache: 



Deutsch 



Deutsch 



(30) Angaben zur Prioritat: 

10 2004 002 929.6 14. Januar 2004 (14.01.2004) DE 

(71) Anmelder (fiir alle Bestimmungsstaaten mit Ausnahme 
von US): DEGUDENT GMBH [DE/DE] ; Rodenbacher 
Chaussee 4, 63457 Hanau (DE). 



(72) Erfinder; und 

(75) Erfinder/Anmelder (nur fur US): WENIGER, Kirsten 

[DE/DE]; Bahrenf elder Chaussee 30a, 22761 Hamburg 
(DE). CAPPIUS, Hans- Joachim [DE/DE]; Wiclefstr. 48, 
10551 Berlin (DE). MULLER, Gerhard [DE/DE]; An 
der Rehwiese 8, 14129 Berlin (DE). 

(74) Anwalt: STOFFREGEN, Hans-Herbert; Friedrich- 
Ebert-Anlage lib, 63450 Hanau (DE). 

(81) Bestimmungsstaaten (soweit nicht anders angegeben, fiir 
jede verfugbare nationale Schutzrechtsart): AE, AG, AL, 
AM, AT, AU, AZ, BA, BB, BG, BR, BW, BY, BZ, CA, CH, 
CN, CO, CR, CU, CZ, DE, DK, DM, DZ, EC, EE, EG, ES, 
FI, GB, GD, GE, GH, GM, HR, HU, ID, IL, IN, IS, JP, KE, 
KG, KP, KR, KZ, LC, LK, LR, LS, LT, LU, LV, MA, MD, 

[Forts etzung aufder ndchsten Seite] 



(54) Title: METHOD FOR DETERMINING COLOUR PERCEPTION IN MULTI-LAYER SYSTEMS 

(54) Bezeichnung: VERFAHREN ZUR BESTIMMUNG DER FARBWAHRNEHMUNG BEI MEHRSCHICHTSYSTEMEN 



Ani=-i?frppi«(riklor g 



TT 



Fteniir-ston einci'otJlisch (fichten Prate 
:Fpr PrOtandif-JtH (f 



K orre ktu r von jj, zxi fi__. 
Inverse Mpni*>-C?.ftu-SmmldlSuii 



."ZT 



Phase nluriklKin 
Piobendicksd 0 \ 
Djecli n litigsi ndex n / 
Messrjeomelne 



Gti^ukoeffizpent ;i, 
Aoir-ctropiefckloi g 



Vorwatxs- Wonrs-Cario-simuraiion 



CFhHsenhrnklii... 
MssEgeorriEtrie / 




WlnUffaklorcnanalyse 



T. 



2& 



AA.... DISPERSION COEFFICIENT \i & 
ANISOTROPY FACTOR G 

BB DIFFUSE REFLECTANCE OF AN OPTICALLY 

DENSE SAMPLE OF SAMPLE THICKNESS D p 

CC... CORRECTION OF M A TO p AK 

INVERSE MONTE CARLO SIMULATION 
DC... PHASE FUNCTION 
SAMPLE THICKNESS Dp 

REFRACTIVE INDEX N 
MEASUREMENT GEOMETRY 

EE.... CORRECTED ABSORPTION COEFFICIENT MflK 
DISPERSION COEFFICIENT \i & 
ANISOTROPY FACTOR G 

FF.... FORWARD MONTE CARLO SIMULATION 

GG PHASE FUNCTION 

REFRACTIVE INDEX N 

MEASUREMENT GEOMETRY 

HH.... DESIRED LAYER THICKNESS 

JJ.... DIFFUSE REFLECTANCE 

KK.... SUITABLE CALCULATION ALGORITHMS 

OR MULTI-FACTOR ANALYSIS 

LL.... COLOUR EFFECT 

COLOUR VALUES IN ACCORDANCE WITH SELECTED 
COLOUR SYSTEM 



FARBWIRKUNG 
Failmcra nach jjevrtlbKem Fsitosystem 



(57) Abstract: The invention relates to a method for determining the colour perception of multi-layer dispersive materials or bi- 
^T} ological materials for layer thicknesses that are respectively selective, by the determination of the diffuse reflectance based on the 
tf^ respective intrinsic optical parameter using Monte Carlo simulations and taking into consideration the measurement geometries, 
anisotropy and the dispersion phase function in order to correctly take into account the multiple internal dispersion of the material. 
The colour effect is calculated from the diffuse reflectance based on the various colour systems in accordance with different algo- 
rithms. 
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(57) Zusammenfassung: Die Erfindung betrifft ein Verfahren zur Bestimmung der Farbwahrnehmung von aus mehreren Schichten 
aufgebauten streuenden Werkstoffen oder biologischen Materi alien bei jeweils wahlbaren Schichtdicken mit Ermittlung der Remis- 
sion auf Basis der jeweiligen intrinsischen optischen Parameter mittels Monte-Carlo-Simulationen unter Einbeziehung der Messgeo- 
metrien, der Anisotropic und der Streuphasenfunktion zur korrekten Berucksichtigung der internen Materialvielfachstreuung. Dabei 
wird die Berechnung der Farbwirkung aus der Remission abhangig von den verschiedenen Farbsystemen nach unterschiedlichen 
Algorithmen durchgefiihrt. 
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Beschreibung 

Verfahren zur Bestimmung der Farbwahrnehmung bei Mehrschichtsystemen 

Die Erfindung bezieht sich auf Verfahren zur Bestimmung der Farbwirkung von streu- 
enden Materialien wie Werkstoffen oder biologischen Substanzen eines Mehrschicht- 
systems, insbesondere einer Schichtabfolge in Zahnen oder dentalen Werkstoffen. 

Ein Ziel der Bestimmung der Farbwahrnehmung ist eine Anpassung oder ein Abgleich 
der Farbwirkung von neuen Mehrschichtsystemen an vorhandene Mehrschichtsysteme, 
beispielsweise als MaBnahme der Qualitatssicherung oder bei der Entwicklung und Be- 
urteilung neuer Materialien und deren Kombinationen. Dies gilt fur Mehrschichtsysteme 
im Bereich von kosmetischen oder asthetischen medizinischen Anwendungen, wie z. B. 
fur Zaknersatz, bis hin zu technischen Gebieten, wie z. B. fur Autolacke oder Kunststof- 
fe. Die Farbwirkung eines Materialsystems ist abhangig von der Schichtdicke und der 
Riickstreuung (Remission) des Lichts und kann gemaB DIN 5033 nach dem Gleich- 
heits-, Dreibereichs- oder Spektralverfahren nach unterschiedlichen Farbsystemen, wie 
z. B. CIELAB oder CIELUV, ermittelt werden. Genutzt werden beispielsweise Messge- 
rate wie Farbspektrometer oder Colorimeter sowie eine visuelle Abmusterung bei- 
spielsweise mittels spezieller Farbskalen, wie diese z. B. fur den Bereich Zahnmedizin 
in der DE-A-196 46 923, der DE-A-101 21 553 oder der DE-A-100 23 840 beschrieben 
werden. 

Aufgrund des bei Mehrschichtsystemen durchzufuhrenden technischen Herstellungs- 
prozesses sowie der Kombination von Einzelschichten kann die Farbwirkung nicht di- 
rekt verglichen, sondern erst im nach hinein, also nach Fertigstellung bestimmt und be- 
urteilt werden, wenn keine Farbkorrektur mehr moglich ist. 



BESTATIGUNGSKOPIE 
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Um eine gewunschte Farbwirkung zu erreichen und die jeweilige Farbwirkung bei vari- 
ierenden Schichtdicken zu bestimmen oder vorauszusagen, sind daher meist zeit- und 
kostenaufwandige Mess- und Versuchsreihen seitens des Herstellers erforderlich, auch 
nach dem Prinzip trial und error. Um dies zu vermeiden, wird versucht, die Farbwirkung 
mit Hilfe von Simulationen zu berechnen. Dies kann fur Einzelschichten und weniger 
komplexe Materialsysteme, die aufgrund ihrer als relativ homogen anzusehenden inne- 
ren Struktur nur eine molekulare Streuung aufweisen, mit ausreichender Genauigkeit 
durchgefuhrt werden. 

Eine solche Berechnung zur Vorhersage der Farbwirkung fur unterschiedliche Schicht- 
dicken wird - sowohl auf dem Gebiet der zahnmedizinischen Wissenschaft als auch im 

► technischen Bereich wie in der Farbindustrie als Basis von, beispielsweise konzentrati- 

onsabhangigen, Farbrezeptierungsberechnungen - bisher mit Hilfe der als Stand der 
Technik geltenden und verwendeten Kubelka-Munk-Gleichung als eindimensionalem 
Modell zur Losung der Strahlungstransportgleichung durchgefuhrt. Aus an Materialpro- 
ben gemessenen Transmissions- oder Remissionseigenschaften werden damit der Ab- 
sorptionskoeffizient Arm und Streukoeffizient Skm berechnet. Auf Basis von Akm und 
Skm konnen die Remissionsspektren fur verschiedene Schichtdicken vorhergesagt wer- 
den, beispielsweise vor weiBem oder schwarzem oder vor nach einem Farbsystem defi- 
niertem Hintergrund. Aus den Remissionsspektren lasst sich die Farbwirkung je nach 
angewendetem Farbsystem berechnen. Beispielsweise kann die Berechnung der Farb- 
wirkung nach CIELAB nach der Norm DIN 5033 durch Umrechung der Remissions- 

i werte in X, Y, Z-Normfarbwerte und nachfolgend in CIELAB-Werte ermittelt werden. 

Diese Umrechnung von gemessenen Remissionswerten in CIELAB-Werte ist auch in 
die Software von Farbspektrometern, wie sie in der Farbindustrie angewendet werden, 
integriert. Durch Vergleich der CIELAB-Werte und Berechnung des Farbabstands AE a b 
nach DIN-Norm 6074 konnen Farbunterschiede zwischen Proben der untersuchten Ma- 
terialien ermittelt werden. 

Mit der Kubelka-Munk-Gleichung konnen jedoch die Messgeometrien der fur die 
Transmissions- oder Remissionsmessungen verwendeten Messgerate nicht berucksich- 
tigt werden und es werden vereinfachende, der Genauigkeit abtragliche, Annahmen ge- 
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troffen. Damit entstehen systematische Fehler, die sich in sichtbaren Farbunterschieden 
zwischen Simulationsergebnis und Wirklichkeit bemerkbar machen konnen. Der groBte 
Nachteil der Kubelka-Munk-Theorie besteht jedoch darin, dass nur Vorgange der Ober- 
flachenre flexion beriicksichtigt werden konnen. Dies wird in Fig. 1 beispielhaft am 
Lichtweg eines Photons 1 nach Reflexion des Lichts an der Oberflache einer Material- 
schicht 2 vorgegebener Dicke dargestellt ist. 

Mit der Kubelka-Munk-Gleichung konnen bei Mehrschichtsystemen, bei denen nicht 
molekulare, sondern strukturelle Streugeometrien vorliegen, die tatsachlich vorliegen- 
den Verhaltnisse nicht mehr korrekt wiedergegeben werden. Denn hierbei kommt es 
durch den Schichtaufbau nnd die Strukturbesonderheiten zu einer Lichtausbreitung mit 
internen Vielfachstreuvorgangen im Material oder Werkstoff. Die Remission und damit 
die Farbwirkung ergeben sich aus einem sehr komplexen Zusammenspiel der optischen 
Eigenschaften der unterschiedlichen Inhaltsstoffe und Komponenten der Schichten auf- 
grund der Lichtausbreitung durch Streuung, Absorption und Brechung an Grenzschich- 
ten. Fig. 2 zeigt beispielhaft den Lichtweg eines Photons 1 sowie einen alternativen 
Lichtweg 1 ' nach multiplen Streuereignissen innerhalb mehrerer hintereinander liegen- 
der Materialschichten 2, 3, 4 unterschiedlicher Dicke und Art. Bei Anwendung der Ku- 
belka-Munk-Gleichung sind daher fur Mehrschichtsysteme mit Strukturstreuung we- 
sentliche Abstriche bei der Genauigkeit der Vorhersage der Remission und somit der 
Farbwirkung hinzunehmen, die sich als visuell wahrnehmbare Unterschiede bemerkbar 
machen. Zur Verbesserung der Vorhersagegenauigkeit der Remission und damit der 
Farbwirkung ist daher die Entwicklung eines Verfahrens erforderlich, das diese struktu- 
rell bedingte Vielfachstreuung korrekt und genau beriicksichtigen kann. 

Der vorliegenden Erfindung liegt das Problem zu Grunde, dass die Farbwahrnehmung 
fur Mehrschichtsysteme aus Kombinationen verschiedener streuender Werkstoffe oder 
biologischer Materialien aus zusammengesetzten, unterschiedlichen Schichten mit ver- 
schiedenen optischen Eigenschaften fur variierende Schichtdicken genauer als bisher 
berechnet und damit vorhergesagt werden kann, ohne jeweils immer wieder Proben, 
bestehend aus den kombinierten interessierenden Schichtdicken, herstellen und deren 
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Farbwirkung, beispielsweise in konventionellen Farbspektralphotometern, messen zu 
miissen. 

Insbesondere soli die Farbwahrnehmung fur Mehrschichtsysteme im dentalen Bereich 
berechnet bzw. vorhergesagt werden, wobei Schichtabfolgen in Zahnen wie Scbmelz 
und Dentin in beliebiger Schichtstarke und dentale Werkstoffe wie Komposits und Ke- 
ramiken von besonderem Interesse sind. 

Zur Losung des Problems sieht die Erfindung im Wesentlichen vor, dass mittels Monte- 
Carlo-Simulationen auf der Basis jeweiliger intrinsischer optischer Parameter der unter- 
schiedlichen Materialien des Schichtsystems umfassend Streukoeffizient |i s , Anisotro- 
piefaktor g und korrigierter Absorptionskoeffizient fi a k und unter Einbeziehung von je- 
weils Brechungsindex n, Dicke d der jeweiligen Schicht sowie Streuphasenfunktion der 
einzelnen Materialien die Remission des Mehrschichtsystems berechnet wird. 

Dabei konnen die intrinsischen optischen Parameter auf Basis von spektrometrischen 
Messungen ermittelt oder aus einer Datenbank entnommen werden. 

ErfindungsgemaB wird ein Verfahren vorgeschlagen, mit dem die visuelle Wahrneh- 
mung der Farbe von Mehrschichtsystemen oder die Ergebnisse deren messtechnischer 
Erfassung, wie beispielsweise durch Farbspektrometer, ausreichend genau bestimmt 
bzw. vorhergesagt werden konnen. Es hat sich uberraschenderweise gezeigt, dass bei 
Vorliegen von strukturellen Streugeometrien in Mehrschichtsystemen mit Simulations- 
rechnungen nach dem Prinzip der Monte-Carlo-Simulation als dreidimensionalem L6- 
sungsansatz der Strahlungstransportgleichung und eines anschlieBenden speziellen Kor- 
rekturverfahrens eine korrekte und wesentlich genauere Bestimmung der Farbwahrneh- 
mung als bisher moglich ist. 

Als Ausgangsdaten fur die erfindungsgemaBe Losung dienen die optischen Eigenschaf- 
ten der jeweiligen im Mehrschichtsystem enthaltenen Werkstoffe oder biologischen 
Materialien in Form der intrinsischen optischen Parameter Absorptionskoeffizient 
Streukoeffizient |i s und Anisotropiefaktor g. 
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Zur Ermittlung der intrinsischen optischen Parameter - auch bezeichnet als mikroskopi- 
sche Parameter die unabhangig von der Materialdicke sind, konnen Verfahren zur 
Anwendung gelangen, die aus der Gewebeoptik fur Dosimetrieberechnungen bei medi- 
zinischen Lasertherapien bekannt sind. So kann mit Hilfe einer inversen Monte-Carlo- 
Simulation aus den in einem Ulbricht-Kugel-Spektrometer an Proben mit einer zuvor 
ermittelten geeigneten Probendicke d gemessenen makroskopischen optischen Parame- 
tern wie der diffusen Remission R<i, der totalen Transmission T t sowie der diffusen 
Transmission T d oder der kollimierten Transmission T c eine Berechnung erfolgen. 

Insbesondere ist vorgesehen, dass die intrinsischen optischen Parameter Streukoeffizient 
|i s? (unkorrigierter) Absorptionskoeffizient ^ a sowie Anisotropiefaktor g eines Materials 
auf der Basis der makroskopischen optischen Parameter des Materials in Form von ins- 
besondere diffuser Remission R d , diffuser Transmission T d und/oder totaler Transmissi- 
on T t und/oder kollimierter Transmission T c unter Benicksichtigung der Streuphasen- 
funktion des Materials, Dicke d einer bei der Bestimmung der makroskopischen Para- 
meter verwendeten Schicht des Materials und Brechungsindex n des Materials mittels 
inverser Monte-Carlo-Simulation berechnet werden. 

Sodann kann fur jedes Material auf der Basis der beiden intrinsischen optischen Para- 
meter Streukoeffizient fi s und Anisotropiefaktor g sowie der Remission einer aus dem 
Material bestehenden optisch dichten Schicht einer Dicke do und unter Beriicksichti- 
gung zumindest der Dicke do, der Streuphasenfunktion und des Brechungsindex n des 
Materials mittels inverser Monte-Carlo-Simulation der korrigierte Absorptionskoeffi- 
zient |u a k berechnet werden. Der Streukoeffizient |u s sowie der Anisotropiefaktor g sind 
bereits aus zuvor erfolgter Berechnung bekannt. 

Mit anderen Worten erfolgt erfindungsgemaB eine Korrektur des berechneten Absorpti- 
onskoeffizients ja a mit Hilfe der Remissonswerte einer optisch dichteren Probe der Pro- 
bendicke do des zu untersuchenden Materials zusammen mit den ermittelten Werten 
von ]x s und g in einer erneuten inversen Monte-Carlo-Simulation zum korrigierten Ab- 
sorptionskoeffizient |a a k. Durch diesen Schritt wird die Genauigkeit der Berechnung der 
Farbwirkung wesentlich erhoht. 
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Die Absorption kann umso genauer bestimmt werden, je groBer die Dicke einer Schicht 
ist. Fur die Bestimmung der Ausgangsdaten |i s und g kann jedoch bei der spektro- 
metrischen Messung keine optisch dichte Probe verwendet werden, da fur diesen Simu- 
lationsvorgang ein gewisser Prozentsatz an kollimierter Transmission, also Lichtdurch- 
lassigkeit, erforderlich ist. 

Die Daten zu den intrinsischen optischen Parametern konnen in Erweiterung des Erfin- 
dungsgedankens materialabhangig in einer Datenbank bereitgestellt werden. Mit dem 
Aufbau einer solchen Datenbank, die die intrinsischen optischen Parameter |n a k, \x s und g 
fur unterschiedliche Materialien oder biologische Substanzen oder dentale Materialien 
enthalt, kann auf diese Werte zuruckgegriffen werden, so dass bei Verwendung dieser 
Materialien in einem Mehrschichtsystem fur die Bestimmung von dessen Farbwirkung 
keine spektrometrischen Messungen mehr erforderlich sind. 

Nachfolgend lassen sich auf der Basis von |i a k, \i s und g durch eine Vorwarts-Monte- 
Carlo-Simulation die Remission fur ein Mehrschichtsystem bei einer wahlbaren Wellen- 
lange und jeweils wahlbaren Dicken der im Mehrschichtsystem enthaltenen Schichten 
berechnen. Anstelle der Remission ist auch die Berechnung der Transmission moglich. 
Die wellenlangenabhangigen Werte der Remission konnen mittels geeigneter Algorith- 
men in Farbwerte wie z. B. die Koordinaten nach CIELAB oder anderer Farbsysteme 
transformiert werden. Moglich ist auch die Anwendung der Multifaktorenanalyse zur 
Ermittlung der Farbwirkung nach neuen, alternativen Farbsystemen. 

Die Berechnung der Farbwirkung ist aufierdem fur unterschiedliche und frei wahlbare 
Messgeometrien moglich. 

Mit dem beschriebenen Vorgehen kann damit die Farbwirkung fur ein Mehrschichtsys- 
tem als Kombination von verschiedenen Materialien wie biologische Materialien oder 
insbesondere zahnmedizinische oder -technische Schichtabfolgen fur variierende 
Schichtdicken der Einzel schichten mit sehr viel hoherer Genauigkeit als bisher be- 
stimmt werden, da zum einen die in den unterschiedlichen Schichten durch die interne 
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Strukturstreuung hervorgerufenen multiplen Streuereignisse korrekt berechnet werden 
konnen und die Genauigkeit zum anderen durch die Bestimmung des korrigierten Ab- 
sorptionskoeffizienten (i a k weiter erhoht wird. Somit ist es moglich, die Farbwirkung fur 
Mehrschichtsysteme so genau vorherzubestimmen, dass keine visuellen Unterschiede 
mehr wahrgenommen werden. Damit kann bei bekanntem Materialverhalten schon vor 
der Herstellung eines Mehrschichtsystems dessen Farbwirkung fur wahlbare Schichtdi- 
ckenzusammensetzungen prazise vorhergesagt werden. AuBerdem kann die Messgeo- 
metrie beriicksichtigt werden und somit auch das Ergebnis von schichtdickenabhangi- 
gen Messungen in unterschiedlichen Spektrometern, beispielsweise von Remissions- 
messungen in Farbspektrometern vorhergesagt werden. 

Besonders geeignet ist das Verfahren fur die Anwendung bei Mehrschichtsystemen in 
der Zahnmedizin fur die Ermittlung der Farbe von dentalen Restaurationsmaterialien 
und Zahnfarbe, bei denen aufgrund der zu unterscheidenden sehr feinen Farbabstufun- 
gen von beispielsweise rotlichem, gelblichem oder graulichem WeiB eine besonders 
hohe Genauigkeit erforderlich ist. Sowohl Zahne, als auch Zahnersatzmaterialien zeigen 
einen Schichtaufbau und Komplexe innere Streugeometrien aufgrund ihrer strukturellen 
Zusanimensetzung. Zahne bestehen aus den Zahnhartsubstanzen Schmelz und Dentin 
mit Prismen- bzw. rohrchenartiger Tubulistruktur, Zahnersatzmaterialien werden zur 
Nachahmung des resultierenden optischen Eindrucks iiir eine optimale Asthetik eben- 
falls aus Schichten unterschiedlicher Opazitat zusammengesetzt mit einem komplexen 
Aufbau aus in einer Matrix eingelagerten Fiillkorpern unterschiedlichster Art und For- 
men. 

Das beschriebene Verfahren, das sich dadurch auszeichnet, dass die Berechnung der 
Remission fur Mehrschichtsysteme mit hintereinander liegenden unterschiedlichen Ma- 
terial schichten wahlbarer Schichtdicke und Eigenschaften mittels Monte-Carlo- 
Simulationen auf Basis der intrinsischen optischen Parameter mit einem korrigierten 
Absorptionskoeffizient |i ak? Streukoeffizient |i s und Anisotropiefaktor g sowie des Bre- 
chungsindex n und der Einbeziehung der Streuphasenfunktion zur Beriicksichtigung der 
internen Vielfachmaterialstreuung erfolgt, ist jedoch auf alle anderen Mehrschichtsys- 
teme mit Strukturstreuung, sowohl im Bereich biologischer als auch technischer oder 
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anderer Systeme anwendbar. Dabei kann auch die Farbwirkung von z. B. nicht in einer 
geraden sondem in einer gekxummten Ebene verlaufenden Schichten berechnet bzw. 
vorhergesagt werden. 

Weitere Einzelheiten, Vorteile und Merkmale der Erfindung ergeben sich aus den An- 
spriichen, den diesen zu entnehmenden Merkmalen -fur sich und/oder in Kombination - 9 
sondern auch aus der nachfolgenden Beschreibung von der Zeichnung zu entnehmenden 
Ausfuhrungsbeispielen. 
Es zeigen: 



Fig. 1 eine schematische Darstellung eines Lichtwegs nach Oberflachenreflexi- 

on an einer Materialschicht, 

Fig. 2 schematische Darstellungen von Lichtwegen in Materialschichten, 

Fig. 3 eine Prinzipdarstellung eines Probenraums mit Ulbricht-Kugel, 

Fig. 4 ein Ablaufdiagramm zur Berechnung intrinsischer optischer Parameter, 

Fig. 5 ein Ablaufdiagramm zur inversen Monte-Carlo-Simulation und 

Fig. 6 ein Ablaufdiagramm zur Berechnung der Remission und nachfolgend der 

Farbwirkung eines Mehrschichtsystems. 



Um die Farbwirkung von streuenden Schichten unterschiedlicher Materialien auf der 
Basis von Remissionsspektren zu bestimmen, werden erfindungsgemafl intrinsische 
(mikroskopische) optische Parameter Absorptionskoeffizient (i a? Streukoeffizient ji S5 
Anisotropiefaktor g der Materialien ermittelt - sei es aufgrund von Berechnungen auf 
der Basis von experimentell ermittelten makroskopischen optischen Parametern, sei es 
unter Berucksichtigung von von Datenbanken zu entnehmenden Werten auf deren 
Basis Monte-Carlo-Simulationen durchgefuhrt werden. Bei der Monte-Carlo-Simulation 
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handelt es sich urn ein fur die Losung von Gleichungssystemen und statistischen Pro- 
zessen angewendetes numerisches Verfahren mit hoher Genauigkeit. 

So konnen mittels einer inversen Monte-Carlo-Simulation mit Hilfe einer z. B. von 
Roggan entwickelten Simulationssoftware WinFit 32, Version 5.2. intrinsische optische 
Parameter auf der Basis gemessener makroskopischer optischer Parameter berechnet 
werden. Zur Messung der makroskopischen optischen Parameter kann ein Ulbricht- 
Kugel-Spektrometer in Form eines Zweistrahl-Spektrometers benutzt werden, wie dies 
prinzipiell der Fig. 3 zu entnehmen ist. 

So ist in Fig. 3 ein Probenraum 10 dargestellt, in dem eine Ulbricht-Kugel 12 angeord- 
net ist. Ferner sind Referenzstrahl 14 und Probenstrahl 16 eingezeichnet. Der Referenz- 
strahl 14 wird liber Spiegel S3', S4' und S5 ? in die Ulbricht-Kugel 12 zur Kalibrierung 
auf einen Referenzstandard 18 gefuhrt. Ulbricht-Kugel 12 und Referenzstandard 18 
konnen aus Spektralon bestehen, wodurch eine 99,8 % ige Reflexion des Lichts fiir eine 
nahezu verlustfreie Detektion der Photonen ermoglicht wird. 

Der Probenstrahl 16 gelangt uber Spiegel ST und S2' nach Durchgang durch eine 
wahlbare strahlenfokussierende optische Linse 20 auf eine nicht dargestellte Probe einer 
Dicke d und in die Ulbricht-Kugel 12. Durch die automatisch bei jedem Messwert nach- 
einander folgende Abgabe von Probenstrahl 16 und Referenzstrahl 14 wird eine Nor- 
mierung der Messwerte durchgefuhrt. Daraus resultiert die Bezeichnung Zweistrahl- 
Spektrometer. 

Die Probe wird im Probenhalter entweder am Transmissionsport 22 oder am Remissi- 
onsport 24 der Ulbricht-Kugel 12 befestigt, je nachdem, ob Vorwarts- oder Rxickwarts- 
Streuung gemessen werden soli. Die jeweils gewahlte optische Linse 20 regelt die Fo- 
kussierung des Lichtspots auf Transmissions- und Remissionsport 22, 24. 

Mit einer entsprechenden Anordnung konnen sodann die makroskopischen optischen 
Parameter diffuse Remission R d , diffuse Transmission Td und totale Transmission T t 
gemessen werden. Alternativ kann bei Verwendung eines anderen entsprechend aufge- 
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bauten Ulbricht-Kugel-Spektrometers - oder einem experimentellen Aufbau mit Ul- 
bricht-Kugel - neben R d und T t statt der diffusen Transmission Td auch die kollimierte 
Transmission T c gemessen werden. Die entsprechenden makroskopischen optischen 
Parameter sind wie folgt definiert. Beim Auftreffen eines Lichtstrahls bildet sich eine 
fur das jeweilige Medium spezifische Verteilung der Photonen im Bestrahlungsvolumen 
aus. Ein Teil der Strahlung wird absorbiert. Ein weiterer Teil der Photonen wird im Me- 
dium gestreut. Diese Streuung kann unterteilt werden in Ruckstreuung als diffuse Re- 
mission Rd und in Vorwartsstreuung als diffuse Transmission Td. Bei geringer Schicht- 
dicke werden auch Photonen ohne vorherige Streuereignisse als kollimierte Transmissi- 
on T c transmittiert. Die Summe aus kollimierter und diffuser Transmission ergibt die 
totale Transmission T t (Td + T c = T t ). Diese makroskopischen Parameter sind wellen- 
langen-, material- und schichtdickenabhangig. 

Fur Materialien, die in Schichtsystemen zum Einsatz gelangen, werden entsprechende 
makroskopische optische Parameter bestimmt. Aus den makroskopischen optischen 
Parametern diffuser Remission Rd, diffuser Transmission Td und/oder totaler Transmis- 
sion T t und/oder kollimierter Transmission T c - zwei von den Transmissionsparametern 
mussen bekannt sein (T d + T c = T t ) - werden sodann durch inverse Monte-Carlo- 
Simulation die intrinsischen optischen Parameter des Materials bestimmt, und zwar Ab- 
sorptionskoeffizient ^i a , Streukoeffizient |i s und Anisotropiefaktor g. Dabei werden bei 
der inversen Monte-Carlo-Simulation zusatzlich Streuphasenfunktion des Materials, 
Probendicke d des Materials, Brechungsindex n des Materials und Messgeometrie be- 
riicksichtigt. 

Messgeometrie bedeutet bei der Verwendung eines Ulbricht-Kugel-Spektrometers z. B. 
Probengeometrie, Blendendurchmesser, Kugelparameter, Strahlendivergenz und Durch- 
messer des Lichtspots. Mit anderen Worten werden die bei der Messung vorhandenen 
geometrischen und optischen Verhaltnisse in der Simulation berucksichtigt, so dass sys- 
tematisch Fehler vermieden werden und die Genauigkeit der Simulation erhoht wird. 
Insbesondere konnen durch die Berucksichtigung der Messgeometrie z. B. auch seitli- 
che Verluste von Photonen an Probengrenzflachen mit einbezogen werden. 
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Eine eindeutige Bestimmung der intrinsischen optischen Parameter erfordert die Mes- 
sung von drei unabhangigen Werten wie R d , T t und T c bzw. T d . 

Als Streuphasenfunktion kann die Henyey-Greenstein-Phasenfunktion oder die Rey- 
nolds-McCormick-Phasenfunktion verwendet werden. Andere geeignete Streuphasen- 
funktionen, die fur das jeweilige Material typisch sind, konnen gleichfalls in Betracht 
gezogen werden. 

Die Streuphasenfunktion kann neben einer mathematischen Formel auch aus einer „look 
up table" bestehen, also Bestimmung der Streuwahrscheinlichkeit tiber den Winkel zwi- 
schen der Richtung des Photons zum Streuereignis und der Richtung nach dem Streuer- 
eignis. Damit kann fur jedes Material eine passende Streuphasenfunktion bestimmt 
werden. Dabei kann eine Streuphasenfunktion durchaus fur mehrere Materialien passen, 
also nicht eindeutig auf ein Material beschrankt sein. Die Heney Greenstein-Funktion 
wird beispielsweise gerne fur eine Vielzahl biologischer Gewebe verwendet. Wenn g 
besonders hoch ist, kann es jedoch von Vorteil sein, die Reynolds-McCormick-Funktion 
zu verwenden. 

Die inverse Monte-Carlo-Simulation ist noch einmal der Fig. 5 zu entnehmen, Prinzi- 
piell arbeitet die inverse Monte-Carlo-Simulation derart, dass von intrinsischen opti- 
schen Parameterwerten ausgegangen wird, die sodann auf die gemessenen oder zur Ver- 
fiigung gestellten makroskopischen optischen Parameter zuruckgerechnet werden. Er- 
gibt sich sodann eine Differenz kleiner als eine Fehlergrenze, so werden die zu Grunde 
gelegten intrinsischen optischen Parameter akzeptiert. 

Entsprechend konnen fur in Frage kommenden Materialien die intrinsischen optischen 
Parameter bestimmt werden. 

Soli zum Beispiel die Remission und damit die Farbwirkung eines drei unterschiedliche 
Materialschichten 2, 3, 4 umfassenden Systems gemafi der Fig. 2 bestimmt, d. h. vor- 
ausgesagt werden, werden die intrinsischen optischen Parameter eines jeden Materials 
sodann einer Vorwarts-Monte-Carlo-Simulation zugrunde gelegt. Ferner werden bei 
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dieser Simulation die Messgeometrien sowie fur jedes Material jeweils die Streuphasen- 
funktion der Materialien und der Brechungsindex n der Materialien beriicksichtigt. Fer- 
ner gehen die Dicken d], d2, d3 der Schichten 2, 3, 4 ein. Die Dicke einer jeden Schicht 
an sich ist dabei konstant. 

Aus der Vorwarts-Monte-Carlo-Simulation ergibt sich sodann eine wellenlangenabhan- 
gige Remission, aus der durch geeignete Rechenalgorithmen oder Multifaktorenanalyse 
die Farbwirkung, also Farbwerte nach gewahltem Farbsystem wie CIELAB berechnet 
werden. 

Um eine hohere Genauigkeit zu erzielen, ist erfindungsgemaB vorgesehen, dass der Ab- 
sorptionskoeffizient korrigiert wird; denn bei den gemaB Fig. 3 durch gefuhrten Mes- 
sungen miissen die Proben eine Dicke aufweisen, die eine Transmission von Licht er- 
moglicht. Die Absorption kann jedoch umso genauer bestimmt werden, je groBer die 
Dicke einer Schicht ist. Daher werden der Streukoeffizient |n s und der Anisotropiefak- 
tor g, die mittels inverser Monte-Carlo-Simulation entsprechend der Fig. 4 und 5 be- 
rechnet worden sind, sowie die Remission einer optisch dichten Probe der Probendicke 
do einer weiteren inversen Monte-Carlo-Simulation zugrunde gelegt, aufgrund der nun- 
mehr ein korrigierter Absorptionskoeffizient ermittelt wird. Bei der Simulation wer- 
den ebenfalls Streuphasenfunktionen des Materials, Probendicke do, Brechungsindex n 
des Materials sowie die Messgeometrie beriicksichtigt. Dies ergibt sich aus dem oberen 
Teil der Fig. 6. 

Fur ein Mehrschichtsystem wird infolge dessen die im Bereich 26 des Ablaufdiagramms 
gemaB Fig. 6 prinzipiell dargestellte Berechnung entsprechend der Anzahl der Schich- 
ten bzw. unterschiedlichen Materialien durchgefuhrt, um sodann gemaB zuvor erfolgter 
Erlauterungen auf der Basis der entsprechenden materialspezifischen intrinsischen opti- 
schen Parameter, also korrigierter Absorptionskoeffizient |a ak , Streukoeffizient |i s und 
Anisotropiefaktor g der Materialien die Remission des Schichtsystems durch Vorwarts- 
Monte-Carlo-Simulation zu berechnen, wobei von den verschiedenen Materialien deren 
Streuphasenfunktion, Brechungsindizes n x sowie Schichtdicken d x und Anzahl x der 
Schichten beriicksichtigt werden, gegebenenfalls erganzend die Messgeometrie. 
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Somit wird die Remission des Gesamtsystems aller Schichten in ihrer geometrischen 
Abfolge berechnet, so dass erwahntermaBen die Phasenfunktionen und Brechungsindi- 
zes der Materialien der Schichten gewunschter Schichtdicken und Anzahl der Schichten 
zu berucksichtigen sind. Die Farbwirkung ergibt sich sodann aus der Remission nach z. 
B. DIN 5033 Farbmessung Teil 1 - 9, Beuth Verlag Berlin, 1979 - 1992. 
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Anspruche 

1. Verfahren zur Bestimmung der Farbwirkung von streuenden Materialien wie 
Werkstoffen oder biologischen Substanzen eines Mehrschichtsystems, insbeson- 
dere einer Schichtabfolge in Zahnen oder dentalen Werkstoffen, 
dadurch gekennzeichnet, 

dass mittels Monte-Carlo-Simulationen auf der Basis jeweiliger intrinsischer op- 
tischer Parameter der unterschiedlichen Materialien des Schichtsystems umfas- 
send Streukoeffizient |i s , Anisotropiefaktor g und korrigierter Absorptionskoef- 
fizient jn a k und unter Einbeziehung von jeweils Brechungsindex n, Dicke d der 
jeweiligen Schicht sowie Streuphasenfunktion der einzelnen Materialien die 
Remission des Mehrschichtsystems berechnet wird. 

2. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

dass die intrinsischen optischen Parameter auf Basis von spektrometrischen 
Messungen ermittelt werden. 

3. Verfahren nach Anspruch 1, 
dadurch gekennzeichnet, 

dass die intrinsischen optischen Parameter aus einer Datenbank entnommen 
werden. 
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4. Verfahren nach zumindest Anspruch 1 , 
dadurch gekennzeichnet, 

dass die intrinsischen optischen Parameter Streukoeffizient fj s , unkorrigierter 
Absorptionskoeffizient \i a sowie Anisotropiefaktor g eines Materials auf der Ba- 
sis der makroskopischen optischen Parameter des Materials in Form von diffuser 
Remission Rd sowie diffuser Transmission T d und/oder totaler Transmission T t 
und/oder kollimierter Transmission T c unter Berucksichtigung der Streuphasen- 
funktion des Materials, Dicke d einer bei der Bestimmung der makroskopischen 
Parameter verwendeten Schicht des Materials und Brechungsindex n des Materi- 
als mittels inverser Monte-Carlo-Simulation berechnet werden. 



5. Verfahren nach zumindest einem der vorhergehenden Anspruche, 
dadurch gekennzeichnet , 

dass fur jedes Material auf der Basis der intrinsischen optischen Parameter 
Streukoeffizient |n s > unkorrigierter Absorptionskoeffizient |i a , Anisotropiefaktor 
g sowie Remission einer aus dem Material bestehenden optisch dichten Schicht 
einer Dicke d D und unter Berucksichtigung zumindest der Dicke do, der Streu- 
phasenfunktion und des Brechungsindex n des Materials mittels inverser Monte- 
Carlo-Simulation der korrigierte Absorptionskoeffizient |i ak berechnet wird. 

6. Verfahren nach zumindest einem der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, 

dass fur das aus unterschiedlichen Materialien bestehende Schichtsystem auf der 
Basis des korrigierten Absorptionskoefffizienten des Streukoffizienten |i s 
und des Anisotropiefaktors g eines jeden Materials unter Berucksichtigung von 
zumindest der Streuphasenfunktion, des Brechungsindex n und Dicke d jeder 
Schicht sowie Schichtreihenfolge mittels Vorwarts-Monte-Carlo-Simulation die 
Remission des Schichtsystems berechnet wird. 
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7. Verfahren nach zumindest einem der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, 

dass bei der Berechnung der intrinsischen optischen Parameter mittels der Mon- 
te-Carlo-Simulation bei der experimentellen Bestimmung der makroskopischen 
optischen Parameter Messparameter und/oder Messgeometrien beriicksichtigt 
werden. 

8. Verfahren nach zumindest einem der vorhergehenden Anspriiche, 
dadurch gekennzeichnet, 

dass aus der Remission die Farbwirkung berechnet wird. 

9. Verfahren nach zumindest einem der vorhergehenden Anspruche, 
dadurch gekennzeichnet, 

dass die Berechnung der Farbwirkung aus der Remission mittels Algorithmen 
oder Multifaktorenanalyse erfolgt. 

1 0. Verfahren nach zumindest einem der vorhergehenden Anspruche, 
dadurch gekennzeichnet, 

dass die Berechnung der Farbwirkung unter Berucksichtigung geometrischer 
Erstreckung wie Krummung des Schichtsystems erfolgt. 
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